---^ Recent results from Pb-Au collisions at 158 GeV/c per nucleon 
G^ ; obtained with the CERES spectrometer 
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j^ ■ During the 1996 lead run time, CERES has accumulated 42 million events, corre- 

sponding to a factor of 5 more statistics than in 1995 and 2.5 million events of a special 
photon-run. We report on the results of the low-mass e"'"e~-pair analysis. Since the most 
critical item is the poor signal-to-background ratio we also discuss the understanding of 
this background, in absolute terms, with the help of a detailed Monte Carlo simulation. 
We show preliminary results of the photon analysis and summarize the results of the 
hadron analysis preliminarily reported on already at QM'97 [0]. 

1. Introduction 

Since electromagnetic probes have a mean free path which is much larger than the size 
of the reaction volume of a heavy ion collision, they carry direct information also from the 
earliest stages of the collision. The CERES spectrometer is optimized to measure low- 
mass e"'"e~-pairs, including the low pair-p^-region which was not accessible in previous 
dimuon experiments at CERN. Also photons can be detected by the conversion method 
(7 — s> e"'"e~). The main advantage of dilepton pairs over photons is that they carry in- 
variant mass which helps to distinguish between different production processes, and in 
particular to observe vector mesons by their direct decay. 
In order to extract signals from the early stages of the collision, detailed knowledge of 



the electromagnetic decays of hadrons after freezeout is necessary. Such a reference was 
provided by measuring the inclusive pair production in p-Be collisions [Q] . Together with 
an electromagnetic calorimeter (TAPS), we have exclusively reconstructed the 7r°- and 
?7-Dalitz decays. It was found that the inclusive pair yield per charged particle can be 
fully understood by hadron decays in p-Be collisions as well as in p-Au collisions. It is 
with respect to this reference that CERES has observed an enhanced pair production, 
first in S-Au collisions and later also in Pb-Au collisions |^. The inclusive photon 
analysis of the S-Au data did not show a significant excess compared to the expectations 
from hadron decays 0. 



2. Experimental Setup 

A main experimental problem a dilepton experiment in nucleus-nucleus collisions has to 
deal with, is the small amount of leptons compared to the huge number of charged hadrons 
(Ne+e-/Nc/i ~ 10~^). This problem is attacked by using two Ring Imaging CHerenkov 
detectors with a 74^ = 32, such that 95% of all charged hadrons do not produce Cherenkov 
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Figure 1. Experimental setup of the CERES experiment during the Pb data taking. The 
dashed curves show the axially-symmetric magnetic field lines. 



light. As shown in Fig. |l|, these RICH detectors are assisted by three external tracking 
detectors, two Silicon Drift Detectors (SDD1,2) 10 cm behind the target and a Multi- 
Wire Proportional Chamber (MWPC) at the end of the spectrometer. The signatures 
of charged-particle tracks are shown in Fig. ^. While the 95 % of all charged hadrons 
below the Cherenkov threshold are only visible in the external tracking detectors, the 
pions with a momentum larger than 4.5 GeV/c are also visible in the RICH detectors 
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Figure 2. Signature of different physics sources in the CERES spectrometer. 

Left Side: Hadrons with a momentum below the Cherenkov threshold (labelled h"*") are 

only visible in the external tracking detectors (SDD's and MWPC). High momentum pions 

(labelled vr"*") and electrons (labelled e"*" and e~) are also visible in the 2 RICH's and can 

be distinguished by the ring radius. The two electron tracks show the typical signature of 

a 7r°-Dalitz decay where the matching ring in RICH2 is not found; two incomplete pairs 

like that are the main source of combinatorial background. The genuine high mass signal 

pairs we are interested in have two well sepaprated tracks already in the first detectors 

(not contained in the figure). 

Right side: The signature of a photon in our special photon run is shown. Since here the 

background is negligible, it is not indicated. 



where they produce a ring the radius of which depends on the relativistic 7 factor of the 
pion and is significantly smaller than the asymptotic radius. Therefore, pion rings are 
well distinguished from electron/positron rings, which always have asymptotic radii. 
The large number of photons per e"'"e~-pair, also of order 10^ demanded to strictly mini- 
mize the radiation length within the detector acceptance (X/Xq < 1 %), so that only 1% 
of the photons convert into a e"'"e~-pair. This is still, together with the decay tt" — >■ e"'"e~7, 
the dominant source for dielectron pairs, both producing pairs with a mass smaller than 
typically 200 MeV/c^ (N,+,-(m>200MeV/c2 )/N,+,-(m<200MeV/c2 ) ~ 10"^). Since it 
is impossible to observe on top of such sources an additional source from an earlier stage 
of the collision, we are mainly interested in pairs with mass larger 200 MeV/c^. 
Because of our limited track reconstruction efficiency and our limited acceptance, it occurs 
very often that from a physical pair only one leg is reconstructed (as indicated in Fig. 0). 
Two of such partially reconstructed pairs in the same event produce a combinatorial 
background pair. To avoid this background a more efficient recognition especially of pairs 
with mass smaller 200 MeV/c^ is required. Therefore the third main idea of the setup 
is to have three detectors (2 SDD's and RICHl) which are able to observe the original 



small opening angle of these pairs which is later used for background rejection. Therefore 
the shape of the magnetic field is such that up to RICHl the spectrometer is field-free, 
between RICHl and RICH2 the charged particles get a deflection in azimuthal direction 
and afterwards the field- lines point to the target (see Fig. |I|), such that the particles have 
again straight tracks. The momentum resolution due to detector resolution and multiple 
scattering is Ap/p = ((2.3% ■ p)^ + (3.5%)2 )V2^ with p in GeV/c. 
For a 5 days running period we took data with a slightly changed setup: between the 2 
SDD's we put a 100 /xm steel foil as photon converter. Photons have a very characteristic 
signature in this setup (see Fig. ^: they have no hit in the first SDD, a double minimum 
ionization pulse hight signal in the second SDD, a double ring in the first RICH and 
two well separated rings in the second RICH accompanied by corresponding hits in the 
MWPC. 

3. Contributions from hadron decays 

As already mentioned the contributions of hadron decays to the dilepton and photon 
yield can be well understood in proton-induced collisions. They depend on the yield 
of the hadrons, the decay branching ratio and, since we have a limited kinematic and 
geometrical acceptance, also on the p_L and rapidity distributions of the hadrons. Except 
for the branching ratios all other ingredients are changing when going from proton-induced 
to lead-induced collisions, as observed by several CERN-SPS experiments: 

The production ratios of the heavier hadrons are enhanced as for example the tj/tt^ 
and the 20/(7r"'"-|-7r~) ratio [0,|1. These particle ratios are best described in a 
thermal model, assuming chemical equilibrium 0, and presently unknown ratios 
have been extrapolated from it. 
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• The inverse m^-slope is systematically increasing with the mass of the hadron [IC 
which is a hint for radial flow. The functional behaviour is: 

inverse m^-slope (mass/ja^ron) = 0.175 MeV + 0.115 MeV/GeV • ma.sShadron 

• While the widths of the hadron rapidity distributions decrease with particle mass 
in proton induced collisions, this is not observed in lead induced collisions [|Tl|] . 

All these observed changes are implemented in the generator simulating the contributions 
of hadron decays and all the data shown in this paper are compared to the outcome 
of this generator. Compared to the pure p-p cocktail used previously this generator 
predicts a 30% larger yield for rHe+e- > 0.2 GeV/c^ (see Fig. ^). For the photon yield the 
correct description of the 7r° distributions in our acceptance is most important, since the 
71-0 _^ ^^ decay dominates the yield (see Fig. |). 

4. Low-mass e"'"e~-pair analysis 

The main challenge in the low-mass pair analysis is to reject the partially recognized 
pairs with small opening angle (as indicated in Fig. |^), while keeping a high reconstruction 



efficiency for pairs with large opening angle []12| . To understand the sources of background 



quantitatively, we have done a detailed MC simulation using a generator for the composi- 
tion of the sources and the decay kinematics. The trajectories of these particles are traced 



through a GEANT simulation of our spectrometer. The resulting dE/dx and Cherenkov 
photon signal is then digitized in the same way as the data. The results are overlaid on 
real events, such that they can be analyzed in the same way as the data. Doing this 
for conversions and 7r°-Dalitz decays on one hand and for pairs with masses larger 200 
MeV/c^ on the other hand, one obtains an objective tool to optimize the cuts. 
The rejection power of these cuts can be seen from Fig. 13 where the absolute contributions 
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Figure 3. Relative num- 
ber of tracks from various 
sources of combinatorial 
background shown along 
the sequence of rejection 
steps applied in the anal- 
ysis. The points indi- 
cate the observed back- 
ground in the data, the 
lines show the estimates 
by the Monte Carlo Sim- 
ulation. 
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of all background sources are followed in chronological order over the different rejection 
steps. The most important rejection steps are called "dE/dx" and "conversion" which 
use the property of the two SDD's to show a double amplitude for two very close tracks 
and "dalitz" which tries to identify close pairs by a close partner ring in RICHl. The MC 
simulation is compared with the observed background in the data. The data follow very 
well the trend of the simulation even in absolute terms. 

At the end of all steps the signal-to-background ratio for pairs with mass larger 200 
MeV/c^ is 1/13, with a mean reconstruction efficiency of 13%. The signal is extracted by 
subtracting the number of like-sign pairs from the number of unlike-sign pairs. By doing 
this one has to be sure that the combinatorial background is equal in these two samples. 
This has been verified on the data by reversing all rejection cuts, such that one obtains a 
"clean" background sample. This study indicates, that we can exclude an asymmetry of 
more than 1% with a 90% confidence level. 



4.1. Results 

Fig. ^ shows the resulting invariant mass distribution of the observed pair yield in 
comparison to the expectations from hadron decays. As previously reported by CERES 
we again observe an enhanced dilepton production in the mass window starting around 
200 MeV/c^ up to the p/uj peak. For comparison the earlier results @] are also shown. 
The centrality selection of the two data sets is slightly different, the mean number of 
charged particles per unit of rapidity is <^ch > = 250 and 220 for '96 and '95, respec- 
tively. The enhancement factor defined by ^^r^e'sured/^Cdronic sources within the mass range 
0.25 < m < 0.7 GeV/c^ is 2.6 ± 0.5 (stat.) ± 0.6 (syst.) for the '96 data. 
The pair yield as a function of the accompanying number of charged particles may be used 



Pb-Au 158AGeV 

p^>200MeV/c 
&^ > 35 mrad 
2.1 <ri< 2.65 

(N , > = 250 / 220 




1.4 1.6 1.8 2 
m^^ (GeV/c') 

Figure 4. Inclusive e"'"e~ mass spectrum 
obtained from the '95 and '96 data anal- 
yses. The data are compared to the sum 
of the expected contributions from hadron 
decays (solid hne). 
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as a tool to disentangle different production mechanisms for the pairs. This dependence is 
shown in Fig. § where the ratio between the measured and generator pair yield is shown 
as a function of multiplicity. Since the expectations from decays of produced hadrons 
grow linearly with the multiplicity, a constant enhancement indicates a linear dNch/d?7 
dependence. Such 'normal' behaviour is displayed by the Dalitz decays dominating the 
two lowest mass bins in Fig. ^, where we do not observe an enhancement. In contrast, 
in the mass bin 0.25 MeV/c^ < m < 0.68 MeV/c^ where the strongest enhancement is 
observed, the pair yield grows stronger than linear. 

We also observe that the enhancement is most pronounced for pairs with a small trans- 
verse momentum (see Fig. |^). While for pairs with pj^ < 500 MeV/c we observe an 
enhancement 5 ± 1.5 for the '96 data, we find only a very modest enhancement under the 
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Figure 6. Inclusive e"'"e~ mass spectra obtained from the '95 and '96 data analyses, se- 
parated into samples with low (left) and high pair-p^ (right). The data are compared to 
the expected contributions from hadron decays. The enhancement is most pronounced 
for low pair-p^. 
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Figure 7. Integrated e"'"e~-pair yield 
as a function of the background-to- 
signal ratio. Shown are the results of 
different analyses using different cut- 
tuning strategies, based either partially 
on data (open symbols) or entirely on 
Monte Carlo simulation (full symbols). 



complementary conditions. This enhancement is consistent with the small enhancement 
in the low-mass region reported by the muon experiments 0JTB| which have no acceptance 
for low p_L. 

In comparison with the previously published data the new high-statistics analysis show a 
statistically consistent but systematically somewhat smaller yield. Therefore we studied 
the stability of the result as a function of the background-to-signal ratio for different cut- 



tuning strategies of the two data samples. The results are summarized in Fig.|^. First of 
all, the results of the '95 sample are seen to vary considerably between the two choices of 
cut-tuning strategies, whether one tunes the cuts partially on the data (which was done 
historically) or with the MC (as described above). These variations are much smaller 
in the '96 sample which is simply a fact of the larger statistics. Much more important 
is, that the signal is stable and does not increase with the B/S ratio. This makes very 
unlikely the hypothesis, that the observed enhancement is due to a incorrectly subtracted 
background. 

Another stability test we did was operating our spectrometer with different polarities of 
the magnetic field. Both polarities show the same result within the statistical errors. 



5. Search for direct photons 

Out of 2.5 million events we reconstructed 34 000 photons which converted into electron 
pairs in the converter foil between the two silicon drift detectors. Due to the clear signature 
shown in Fig. 1^, the background is negligible here ||14 . 
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Figure 8. Inclusive photon p^ 
distribution in comparison with 
the expectations from hadron de- 
cays. The data and the gener- 
ator are consistently normalized 
by the number of photons per 
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the CERES acceptance. The 
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tematic uncertainty of the gener- 
ator, while the brackets indicate 
the systematic error of the data 
points. 
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Before obtaining the p_L distribution one has to correct the data for detector inefficiencies 
and momentum resolution. This correction is obtained with the same MC simulation 
used for the low-mass pair analysis. The corrected p_L-distribution is shown in Fig. ^ 
The lower bound of 400 MeV/c is due to the poor efficiency for very low-momentum 
(large-deflection) tracks, while the upper bound is due to limited statistics. In this in- 
terval the data closely follow the expectations from hadron decays. It is the challenge of 
the photon analysis to do the normalization of both yields as accurately as possible and 
to minimize the systematic errors such that one is sensitive to an additional source on 
the per cent level (corresponding to a low-mass pair enhancement on top of the 7r°-Dalitz 
peak). The normalization is done with respect to the number of photons counted per 
71^ in the same acceptance with a p_L > 400 MeV/c. Therefore we need to know the 
vr'^-px-distribution for our centrality and acceptance. This can be at best estimated from 
our measured h~-spectrum (as described below) by correcting it for contributions from 



p, K~ and secondaries from the t] -^ Svr^ decay. In addition we have also investigated 
charged pion and uncharged pion distributions from other SPS experiments; these serve 
to estimate the systematical errors connected with our choice. 

As indicated in Fig. |] the statistical errors of the measurement are negligible. More 
critical are the systematic errors due to the Monte Carlo correction (6.3 %), the uncer- 
tainty in the absolute normalization (6.6 %) and the uncertainty of the conversion length 
(2.6 %). The expectations from hadron decays also contain a systematic error originating 
from the uncertainty in the p_L distribution and the rj/n^ ratio (4 %). 
Within these errors we obtain a non-significant excess for photons with p_L > 400 MeV/c 
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Within the errors this ex- 



cess is constant as a function of multiplicity. Further studies will verify this result with 
an independent data sample taken with a different magnetic field, to reduce thereby the 
systematic uncertainties. 



6. Hadron analysis 

As indicated in Fig. |], we use two independent methods of reconstructing charged 



hadrons [|15|. First we can reconstruct low-momentum hadrons by matching tracks in the 
SDD's to hits in the MWPC and determine the momentum by the deflection between 
them. This method suffers from the large background (S/B = 2:1), which is determined 
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Figure 9. p^-distributions 

for charged hadrons (labelled h~) 
and identified charged pions (la- 
belled (7r++7r-)/2). The charged 
hadron spectrum is shown before 
(full circle) and after subtracting 
the contributions from K~ and 
anti-protons (open circle), the lat- 
ter being in excellent agreement 
with the overlapping part of the 
pion spectra. 



by the event-mixing method. The second method is sensitive to pions with p > 4.5 
GeV/c, since these pions emit Cherenkov light. The Cherenkov ring radius is used for 
the momentum determination. The p^-slopes of both methods show a very good agree- 
ment in the overlapping region (see Fig. |^). In the region 1.5 < px < 3.5 GeV/c of 
the inclusive spectra, the inverse slope is 245 ± 5 MeV and increases by 2.4 % with in- 
creasing centrality, in the range of the most central 35% of the geometric cross section. 
Identified vr"*" and tt~ show the same slope and the ratio tt~ /tt^ = 1.030 ± 0.005 for 
1 < Px < 2.4 GeV/c. Compared to our previous pion analysis done for S-Au collisions 



10 



the inverse slope as a function of m^ does not follow the systematics of the Cronin effect, 
since the slope at small m^ is much higher. This is a hint for radial flow in Pb collisions 
already observed by other experiments |1^ . 



We also studied anisotropic transverse flow by using only the SDD's in terms of a Fourier 
decomposition of the azimuthal charged multiplicity distribution. Our preliminary results 
show an anisotropic flow up to the 4th harmonic, with the higher harmonics having nearly 
the same strength than the first two. 

7. Conclusion and Outlook 

We have corroborated the low-mass pair enhancement previously observed by CERES. 
This enhancement is most pronounced in the mass window from 200 < m < 700 MeV/c^, 
at low p± and in events with a high charged particle multiplicity. In order to draw final 
conclusions on whether these data exhibit a signal for the restoration of chiral symmetry 
[0 better data are necessary. Therefore CERES has upgraded with a TPC installed at 



the downstream end of the present spectrometer [Q to improve the signal-to-background 
ratio and the mass resolution. This TPC was first tested in 1998 |T^. CERES will con- 
tinue data taking in the next two years with higher statistics and at two different beam 
energies (40 and 158 GeV per nucleon). 
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